Magnetic field-assisted manipulation and entanglement of Si spin qubits 
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Architectures of donor-electron based qubits in silicon near an oxide interface are considered 
theoretically. We find that the precondition for reliable logic and read-out operations, namely the 
individual identification of each donor-bound electron near the interface, may be accomplished by 
fine-tuning electric and magnetic fields, both applied perpendicularly to the interface. We argue 
that such magnetic fields may also be valuable in controlling two-qubit entanglement via donor 
electron pairs near the interface. 
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Spin qubits in semiconductors (e.g. GaAs, Si) are 
among the most promising physical systems for the even- 
tual fabrication of a working quantum computer (QC). 
There are two compelling reasons for this perceived im- 
portance of semiconductor spin qubits: (i) electron spin 
has very long coherence times, making quantum er- 
ror correction schemes feasible as a matter of principle; 
(ii) semiconductor structures provide inherently scalable 
solid state architectures, as exemplified by the astonish- 
ing success of the microelectronics technology in increas- 
ing the speed and efficiency of logic and memory opera- 
tions over the last fifty years (i.e. 'Moore's Law'). These 
advantages of semiconductor quantum computation ap- 
ply much more to Si than to GaAs, because the electron 
spin coherence time can be increased indefinitely (up to 
100 ms or even longer in the bulk) in Si through isotopic 
purification 0,0, U whereas in GaAs the electron spin co- 
herence time is restricted only to about 10 us [JJ, |2j, |4| . It 
is thus quite ironic that there has been much more sub- 
stantial experimental progress in electron spin and 
charge qubit manipulation in the III-V semiconductor 
quantum dot systems [5| than in the Si:P Kane computer 
architecture p]. In addition to the long coherence time, 
the Si:P architecture has the highly desirable property of 
microscopically identical qubits which are scalable using 
the Si microelectronic technology. The main reason for 
the slow experimental progress in the Kane architecture 
is the singular lack of qubit-specific quantum control over 
an electron which is localized around a substitutional P 
atom in the bulk in a relatively unknown location. This 
control has turned out to be an impossibly difficult exper- 
imental task in spite of impressive developments in mate- 
rials fabrication and growth in the Si:P architecture using 
both the 'top-down' and the 'bottom- up' techniques [2|. 
It is becoming manifestly clear that new ideas are needed 
in developing quantum control over single qubits in the 
Si:P QC architecture. 

In this Letter we suggest such a new idea, establish- 
ing convincingly that the use of a magnetic field, along 



with an electric field, would enable precise identification, 
manipulation and entanglement of donor qubits in the 
Si:P quantum computer architecture by allowing control 
over the spatial location of the electron as it is pulled 
from its shallow hydrogenic donor state to the Si/Si02 
interface by an electric field. Additionally, the magnetic 
field could be used to control the spatial overlap of the 
electronic wavefunctions in the 2D plane parallel to the 
interface (i.e. similar to a MOSFET geometry), thus en- 
abling external manipulation of the inter-qubit entangle- 
ment through the magnetic field tuning of the exchange 
coupling between neighboring spin qubits ||. 

Fig- H] highlights the basic physical effects explored in 
our theoretical study: Substitutional P donors in Si, sep- 
arated by R, are a distance d from an ideally flat Si/Si02 
(001) interface. Uniform electric F and magnetic B fields 
are both applied along z. When a single P atom is con- 
sidered (R — > oo limit), and in the absence of external 
fields, the active electron is bound to the donor poten- 
tial well (Wd) forming a hydrogenic atom. If a uniform 
electric field is also present, an additional well is formed 
which tends to draw the donor electron toward the inter- 
face. The interface well (Wj) is triangular-shaped along z 
while in the xy 2D plane the confinement is still provided 
by the 'distant' donor Coulomb attraction. The ground- 
state wavefunctions for each well are given in Fig. GJb) . 
This corresponds to a particular value of the applied field, 
the critical field which we denote by F = F c (d) , such that 
the expectation values of the energy for the Wi and Wd 
ground-eigenstates are degenerate. Fig. ^c) shows how 
the system in (b) changes due to a magnetic field ap- 
plied along z. The interface ground state wavefunction 
undergoes the usual 'shrinking' perpendicular to B, and 
the energy at Wj is raised, while no significant changes 
occur at Wd- As a consequence, we find that by prop- 
erly tuning B it is possible to control the electron loca- 
tion along the applied magnetic field direction, partially 
or completely reversing the effect of the electric field. 

We base our quantitative description of this problem 



2 



F,B 



(a) 



Si0 2 Si 



?P 




-30 -20 -10 
z (nm) 



FIG. 1: (Color online) (a) Schematic configuration of donors 
in Si near the interface with an oxide barrier under applied 
electric F and magnetic B uniform fields. In (b) and (c), the 
thick lines show the electronic confining double- well potential: 
The interface well is labeled Wi and the donor well, Wd- The 
dashed lines give the decoupled ground-eigenfunctions in the 
wells, and the thin horizontal lines indicate the expectation 
value of the energy in each well. Results for d = 30 nm and 
F = 13.5kV/cm?s F c (d), are given in (b) for 5 = and in 
(c) for B = 10 T. Parameters defining typical wavefunction 
widths (see text) are indicated in (c). Magnetic field effects 
tend to be stronger for larger d, as the interface state is less 
affected by the strongly confining Coulomb potential at Wd- 



on the single- valley effective^mass approximation, leading 
to the model-Hamiltonian 



H = T + eFz- 
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The magnetic field vector potential, A = B (y, —x, 0) /2, 
is included in the kinetic energy term, T = 

J2r,=x, y ,z ^ I ( 2m »?) (id/&n + eA n/ (frc)) 2 ; where the effec- 
tive masses m x — m y = m± = 0.191m and m z — 
rn 1 1 = 0.916 m, account for the Si conduction band val- 
ley's anisotropy. The electric field defines the second 
term in H, the third and fourth terms describe the donor 
and its image charge potentials, while the last term is the 
electron image potential. The fourth term involves the 
lateral coordinate p = (x,y). The image-related terms 
are proportional to Q — [e.2 — ei)/[ei(e2 + ei)], where 
ei = 11.4 and — 3.8 are the Si and SiC>2 static di- 



electric constants. It is convenient to rewrite the kinetic 
energy term as 



T = T Q + X -u? cP 2 
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where Tq is the kinetic energy for B = and oj c — 
eB/(m±c). The second term allows interpreting the 
effect of the magnetic field as providing an additional 
parabolic potential in the xy plane, which increases the 
kinetic energy and enhances the lateral confinement of 
the electron wave- function. The last term in T gives 
zero or negligible contribution for the donor electron low- 
energy states being investigated, and will be neglected 
here. 

We solve the double-well problem described by the 
model Hamiltonian in Eq. in the basis of the un- 
coupled solutions to the individual wells Wi and Wd, 
obtained variationally from the ansatz [IJ : 



*/0,2) = f a (z) x g p {p) 



(3) 
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where a, (3, a and b are the variational parameters. The 
(z + d) factors guarantee that the wavefunctions are zero 
at the interface (z = —d), as we assume that the insulator 
provides an infinite barrier potential. For d > 6 nm and 
B = 0, we find that a and b coincide with the Kohn- 
Luttinger variational Bohr radii for the isolated impurity 
(d — * oo), where a = 2.365 nm and b = 1.36 nm. 

As an illustration of the uncoupled wells variational so- 
lutions, we present in Fig.^the calculated wavefunctions 
for d = 30nm. The dash-dot lines correspond to f a {z) 
and ^^(O, z), while the dashed lines on the left represent 
gp{p) for p = (x,0) [or equivalently (0, y)]. It is clear 
from this figure that the four variational parameters de- 
fine relevant length scales involved in the problem. In the 
presence of a magnetic field, there is an additional length 
scale (the magnetic length) given by Xb — (eB) 
which defines the typical lateral confinement produced 
by the magnetic field alone, regardless of other poten- 
tials in the problem. The magnetic field effect is signifi- 
cant only for values of B > B c such that Xb ^ c p where 
c p is the lateral confinement length in the absence of the 
field. For the donor potential well, c p ~ a = 2.365 nm, 
so that the field required to appreciable affect the Wr>- 
ground state eigenenergy or wavefunction is ^120 
T! On the other hand, for the Wj well, the confinement 
length at the interface for B = is 1//3 ~ 18.5 nm so 
that a much smaller B[ ~ 2 T is sufficient to affect the 
interface ground state. The magnetic field strengths we 
consider here (up to 10 T) are not large enough to affect a 
or b, but important effects are obtained in lateral confine- 
ment at the interface as shown in Fig. ^ by comparison 
of B = in (b) with B = 10 T in (c). 
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FIG. 2: (Color online) (a) Shift in the energy Ei, associated 
with the interface state, versus magnetic field for three values 
of the donor distance to the barrier (d). The shift in energy 
of the donor state (Ed) is negligible compared to these values 
and is not shown, (b) Expectation value of the z-position for 
the electronic ground state versus magnetic field. For each 
d the electric field is tuned to a value just above the critical 
field: F = F c + 60 V/cm. The electron starts at a position 
close to the interface and moves parallel to the magnetic field 
and against the electric field ending up at the donor well. 
The values of the critical electric field at each distance are 
F c (15.8 nm) = 27.44 kV/cm, F c {22 nm) = 18.85 kV/cm, 
F c (30 nm) = 13.44 kV/cm. (c) Same as (b) for F = F C + 120 
V/cm. 



The expectation value of the energy for each well 
ground state, Ej = (^! j\H\^ j) D ), is given by the hor- 
izontal lines in Fig. 2] The critical field condition F = F c 
corresponds to Ej = Ed and is illustrated in Fig.^b). In 
Fig. njc) we also indicate the value of the B = energies 
by a thin dotted line, and we note that under a 10 T 
magnetic field Ej is raised by 2 meV, while Ed undergoes 
a comparatively negligible shift (by 0.13 meV). Fig. Ufa) 
shows the energy shifts, AE(B) = Ei(B) — £j(0), for 
three particular values of d: As expected, the effect of 
the magnetic field is stronger for larger d. 

Consequences of this shift in energy are shown in 
Figs. |2I b) and (c) where we give the expectation value 
of the z-coordinate of the electronic ground state Wo = 
Ci^i + Cd^d- We start with an electric field slightly 
above F c (d), hence, Cj w 1 and Cd ~ 0. Under an 
increasing magnetic field, the energy shift at the inter- 
face eventually detunes Ed and Ei in such a way that 
Cj s» and Cd ~ 1, i-e. the electron is taken back to 
the donor moving parallel to the magnetic field B and 
against the electric field F. How big a magnetic field is 
needed for this process to be completed depends on d and 
on the value of the static electric field. For those donors 
further from the interface, the passage of the electronic 
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FIG. 3: (Color online) (a) Reduction of the lateral confine- 
ment length of the interface state, 1//3, when a perpendicular 
magnetic field is applied. The effect is stronger for donors 
that are further away from the interface, (b) Overlap S versus 
magnetic field for inter-donor separation R = 30nm (which 
corresponds to a planar density n = 10 11 cm -2 ). 



ground state from the Wj to Wd occurs more abruptly 
as a function of B. The effect of the value of the electric 
field is also relevant, as illustrated in Fig. EJb) and (c). 
In Fig. |2^b), F = F c + 60 V/cm and the passage occurs 
around B ~ 2.2 T, while in Fig.^c), F = F; + 120 V/cm 
and the needed field is B ~ 3.2 T. Therefore, fine tun- 
ing of the electric field is required to observe this phe- 
nomenon at reasonably small magnetic fields. We pro- 
pose to use this result as a means of differentiating donor 
electrons from other charges that may be detected 01 m 
the architecture shown in Fig. ^a) . Tunneling times are 
not significantly affected upon magnetic fields considered 
here and remain a function of d alone 0. 

The other significant effect of the magnetic field on 
the interface states is the transverse shrinkage of the 
wave function. Fig. 0{a) shows the lateral confine- 
ment length 1//3 as a function of the magnetic field. 
For d — 15.8 nm, the wave-function width is reduced 
by a factor of 0.8 when a 10 T magnetic field is ap- 
plied, while for d = 30 nm a stronger reduction fac- 
tor of 0.6 is obtained. The significance of this en- 
hanced confinement may be accessed through the over- 
lap between interface electrons coming from neighbor- 
ing donors separated a distance R. The overlap quan- 
tifies whether we can consider these electrons as sep- 
arate entities or they will form a 2-dimensional elec- 
tron gas (2DEG). One obtains S = {gp{p)\gp{p + 
R)) = exp(-/? 2 i?74)exp[-i?7(16A|j/3 2 )] where we 
note that, due to gauge invariance, gp(x + R,y) — 
TT- 1 (3exp{-P 2 [(x + R) 2 +y 2 }/2}exp{iyR/(A\ 2 B )}. The 
results for the particular inter-donor distance R — 30 nm, 
which corresponds to a planar density n ~ 10 11 cm -2 , are 
shown in Fig.^b). We estimate that an overlap S < 0.1 
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FIG. 4: (Color online) (a) Overlap versus interdonor distance, 
(b) and (c) show the magnetic field required to reduce the 
overlap to 0.1 and 0.2 respectively. R = 20 nm corresponds 
to a planar density n — 2.5 x 10 11 cm -2 and R — 40 nm to 
n = 6.25 x 10 10 cm -2 . 



would guarantee that the electrons do not form a 2DEG. 
Fig.QJa) presents the overlap versus inter-donor distance 
for B = 0. The calculated S is generally very large for the 
experimentally reasonable distances d and R considered 
here. For the donors further from the interface (d = 30 
nm), the condition S < 0.1 requires either the applica- 
tion of a magnetic field (around 6 T for R = 30 nm) or 
a smaller planar density n % 3 x 10 10 cm~ 2 (which cor- 
responds to R > 55nm). The magnetic fields required 
to get a reduction of the overlap to S = 0.1 and 0.2 are 
shown in Fig. Qfb) and (c), respectively, as a function 
of the inter-donor distance. For the donors closer to the 
interface (d = 15.8 nm), the effect of the magnetic field 
is not as dramatic: The lateral confinement provided by 
the donor potential is strong enough to give S = 0.1 for 
B = and R ~ 37 nm [n ~ 7.4 x 10 10 cm" 2 ). Note 
that in the case of neighboring donors placed at differ- 
ent distances from the interface d, S(R) would present 
oscillations due to valley interference effects [llj . 

For qubits defined at the donor sites it is certainly 
more reliable to perform operations involving two-qubit 
entanglement at the more directly accessible interface re- 
gion. In this case, depending on R and d, the confining 
potential provided by the donors may need to be com- 
plemented by additional surface gates 01 and/or static 
magnetic fields. Time-dependent magnetic fields can be 
invoked in switching the exchange gates for two-qubit 
operations ||. 

In summary, we demonstrate here that, in addition to 



the control of donor charges by electric fields, as usual in 
conventional semiconductor devices, relatively moderate 
magnetic fields (up to 10 T) may provide relevant in- 
formation and manipulation capabilities as the building 
blocks of donor-based QC architectures in Si are being 
developed. Uniform magnetic fields may displace donor- 
bound electrons from the interface to the donor nucleus 
region, as illustrated in Fig. This effect, which could 
be monitored via surface charge detectors, allows differ- 
entiating donor electrons from spurious surface or oxide- 
region bound electrons, which would not exhibit this type 
of behavior. For electrons originating from neighboring 
donors and drawn to the surface by an electric field, the 
magnetic field lateral confinement effect may provide iso- 
lation between electrons (see Fig.OJ, allowing smaller in- 
terdonor distances to be exploited (see Fig. 0J, as well 
as additional control in two qubits operations via surface 
exchange gates similar to the spin manipulation in quan- 
tum dots in GaAs. Defects such as dangling bonds at 
the Si/Si02 interface are important sources of spin de- 
coherence, thus spin coherence times near the interface 
are expected to be significantly reduced as compared to 
the bulk. Performing two-qubit (spin-spin) entanglement 
near the SiC>2 interface, as suggested here, requires ex- 
tremely careful interface optimization. We believe that 
the use of an external magnetic field along with the FET 
geometry near a Si-Si02 interface should allow in the 
near future significant experimental progress in the cur- 
rently stalled, but potentially important, donor qubit 
based Kane spin QC architecture. 
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